Nonperturbative Scaling Theory of Free Magnetic Moment Phases in Disordered Metals 
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The crossover between a free magnetic moment phase and a Kondo phase in low dimensional disordered 
metals with dilute magnetic impurities is studied. We perform a finite size scaling analysis of the distribution 
of the Kondo temperature as obtained from a numerical renormalization group calculation of the local magnetic 
susceptibility and from the solution of the self-consistent Nagaoka-Suhl equation. We find a sizable fraction 
of free (unscreened) magnetic moments when the exchange coupling falls below a disorder-dependent critical 
value J c . Our numerical results show that between the free moment phase due to Anderson localization and 
the Kondo screened phase there is a phase where free moments occur due to the appearance of random local 
pseudogaps at the Fermi energy whose width and power scale with the elastic scattering rate 1/r. 



The Kondo problem is of central importance for under- 
standing low-temperature anomalies in low-dimensional dis- 
ordered metals CtiHSSHHJiliEIiillil, such as the 

saturation of the dephasing rate II 1311 and the non-Fermi-liquid 
behavior of certain magnetic alloys (3(3. For a clean metal, 
the screening of a spin- 1/2 magnetic impurity is governed by 
a single energy scale, the Kondo temperature Tr. Thermo- 
dynamic observables and transport properties obey universal 
functions which scale with Tk- Thus, in a metal where non- 
magnetic disorder is also present, two fundamental questions 
naturally arise: (i) Is the Kondo temperature modified by non- 
magnetic disorder? (ii) Is the one^parameter scaling behav- 
ior still valid? It is well known Jj, H H 0, [H that mag- 
netic moments remain unscreened when conduction electrons 
are localized due to disorder. However, in weakly disordered 
two-dimensional electron systems, the localization length is 
macroscopically large, and so is the number iV c of eigenstates 
with a finite amplitude at the position of the magnetic im- 
purity. In this case, one does not expect to find unscreened 
magnetic moments for experimentally relevant values of the 
exchange coupling J. 

Another situation where magnetic moments remain free in 
metals at low temperatures occurs when the density of states 
has a global pseudogap at the Fermi energy Ep, namely, 
u(E) ~ (E - E F ) a , where a > 0]. In clean metals, 
the pseudogap quenches the Kondo screening when J falls 
below a critical value J c (a). So far, only a few values of a 
have been realized experimentally: a = 1 in graphene and in 
(i-wave superconductors and a = 2 in p-wave superconduc- 
tors. In this Letter we examine the quantum phase diagram 
of magnetic moments diluted in two-dimensional disordered 
metals using a modified version of the numerical renormal- 
ization group (NRG) method. We find a free moment phase 
which we attribute to the random occurrence of local pseu- 
dogaps. The existence of free moments is confirmed directly 
with NRG by the Curie-like behavior of the the local mag- 
netic susceptibility at low temperatures. Finite-size scaling is 
performed to demonstrate the robustness of our finding. Fur- 
thermore, the distribution of Kondo temperatures obtained nu- 
merically from NRG is found to agree well with earlier results 



based on the solution of the Nagaoka-Suhl equation J5, 11 |. 

We consider the Kondo Hamiltonian of a magnetic impurity 
diluted in a disordered metal, 
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Here, E n are the eigenenergies of the noninteracting conduc- 
tion electrons in the disordered metal, which we describe by 
the Anderson tight binding model, 
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with band width D, nearest-neighbor hopping amplitude t, 
and random site potentials VI, drawn from a flat box distri- 
bution of width W centered at zero. We consider square lat- 
tices of length L with N = L 2 sites and periodic boundary 
conditions. The exchange coupling matrix elements are then 
given by J nn < = J ip* {r)ip n > (r), with ip n (r) being the ampli- 
tude of the single-electron eigenfunction at the position r of 
the magnetic impurity. Energies are given in units of t. 

Finite-Size Scaling. First, we obtain all eigenenergies and 
eigenfunctions of Hy by using state-of-the-art numerical di- 
agonalization techniques for a large set of realizations of 
the potential V [16]. Within one-loop approximation, the 
Kondo temperature is then obtained from the solution of the 
Nagaoka-Suhl equation (NSE) [17], 
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We assume that E-p is between energy levels, with even num- 
ber of electrons in the ground state of the system. For odd 
number of electrons, the single electron at the Fermi energy 
forms a singlet with the magnetic impurity with binding en- 
ergy J. Solving Eq. (O numerically, one finds strong de- 
viations from the Gaussian behavior even for weak disorder 
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FIG. 1 : (Color online) The distribution function of Kondo tempera- 
tures obtained from the NSE [Eq. ([3}] for electrons in tight-binding 
square lattices with disorder W — 3. Ensembles of 100 samples of 
size L = 30, 5 samples for 40, and one sample for L = 50, 70 were 
used. The exchange coupling was fixed to J/D = 0.35. 
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FIG. 2: (Color online) The distribution of Kondo temperatures ob- 
tained with NRG for lattices of size L = 70. (a) W = 2, 
J/D = 0.3. (b) W = 3, J/D = 0.35 for one realization of the 
disorder potential. For comparison we plot the distribution obtained 
from NSE. Hi gher l oop corrections are accounted for by rescaling 
T K with 0.7 y/j/D. The error bars denote the statistical error. 



I^. Figure Q] shows for W = 3 that a double-peak struc- 
ture persists as the lattice size is increased from 900 to 4,900 
sites. This structure evolves into a power-law divergence in 
the strong-disorder limit |H, JJ.]. In order to verify that these 
features are not an artifact of the one-loop approximation, we 
performed a comparative analysis of the Kondo temperatures 
obtained with the nonperturbative NRG method for different 
system sizes and disorder strengths. 

NRG for Disordered Systems. In order to apply the NRG 
method to disordered systems, we use as a basis set the eigen- 
functions and eigenenergies of the the tight-binding model, 
Eq. (O, as obtained by the numerical exact diagonalization. 
We perform the Wilson logarithmic discretization of the con- 
duction band and replace the interval of eigenenergies be- 
tween (D/2)A- n+1 and (D/2)A~ n , n = 0, 1, 2, by one 
energy level equal to the average of eigenenergies in this inter- 
val. The NRG discretization parameter is set to the smallest 
tractable value, A = 1.5. Next, with the help of the well- 
known tridiagonalization procedure, we numerically map this 
model onto a discrete Wilson chain. Then, we use an iterative 
diagonalization procedure, keeping 1,500 states after each it- 
eration [18]. We calculate the local magnetic susceptibility 
Xi oc (T), which is proportional to the local correlation func- 
tion of the impurity spin. The susceptibility is found to have 
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FIG. 3: The local spin susceptibility as function of temperature cal- 
culated with NRG for exchange coupling J/D = 0.35, disorder am- 
plitude W — 2, and lattice size L = 70. (a) At a site where Tk 
is maximal. Dots: Results obtained with the continuous-time Quan- 
tum Monte Carlo method (CT-QMC) (HIH. (b) At a site where 
the magnetic impurity remains free. The arrows indicate the Kondo 
temperature of the clean system for the particular value of J/D. In- 
sets: The absolute value square of the eigenfunction amplitudes at 
those sites as function of energy E, where E — denotes the Fermi 
energy. The mark indicates the mean level spacing A. Energies and 
temperatures are given in units of t. 



a smooth temperature dependence (see Fig. [3]) which can be 
used to define Tjc through the following criterion: Tx\oc{T) 
crosses over from the free spin value 1/4 to a decay linear 
in T as the temperature is lowered below the Kondo temper- 
ature Tk- For a clean flat band system, Xioc(^) coincides 
with the impurity susceptibility Ximp, which is obtained as 
the difference between the susceptibility of the electronic sys- 
tem with and without the impurity 
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the Kondo temperature as the crossover temperature where 
T K Ximp(T K ) = 0.07 USB. For disordered systems, we 
find that Tximp(T) can strongly deviate from the universal 
scaling curve of the clean system ll9l l22Tl . It turns out, how- 
ever, that this is an artifact of the definition of Xi mp : The sus- 
ceptibility of the conduction electrons fluctuates widely and 
can result in negative values of Ximp 12411 . Therefore, we con- 
sider the local magnetic susceptibility as obtained by differ- 
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FIG. 4: (Color online) Fraction of free moments as function of ex- 
change coupling obtained by two methods: (a) From the condition 
that the NSE has no solution; (b) from temperature dependence of 
local spin susceptibility calculated with NRG for disorder W — 2. 
The corresponding NSE result, and the fraction of free moments in 
a clean sample (step functions) are plotted for comparison. Lattice 
sizes are indicated. The lines are guides to the eye. 



entiating the magnetic moment of the impurity with respect 
to a magnetic field that acts only at the impurity B25I1 . In Fig. 
|2l we show the distribution of Kondo temperatures when Tk 
is obtained from the Wilson condition TkXioc{Tk) = 0.07. 
The data was extracted from a single sample of size L = 70 
using two distinct sets of disorder and exchange coupling am- 
plitudes: W = 2 and J = 0.3 L> (Fig. and W = 3 
and J = 0.35 D (Fig. [2J3). For comparison, we also plot 
the distribution of Tk obtained from the solution of the one- 
loop NSE, where we accounted for the known higher-loop 
correction by rescaling T% with 0.7 ■\/J/D. The agreement 
is remarkable and demonstrates that the double peak struc- 
ture is not an artifact of the one-loop approximation. In Fig. 
Hk, we plot the local impurity spin susceptibility, multiplied 
by temperature T, for the site with maximal Tk for a given 
realization of the disorder. In order to check that the modi- 
fied version of NRG is well suited for disordered systems, we 
also plot in Fig. [3^ results obtained with the continous time 
quantum Monte Carlo (CT-QMC) method (dots) |Q. For 
temperatures close to the Kondo temperature both methods 
agree well. The small deviations seen at larger temperatures 
can be attributed to the fact that the average of the wavefunc- 
tion amplitudes in each Wilson discretization interval results 
in stronger suppression of fluctuations at higher energies. We 
note that the local density of states (LDOS), which is shown 
in the inset, has most of its weight in the lower half of the 
band, close to the Fermi energy. In contrast, in Fig. |3J) we 
show the impurity susceptibility times T for a site where the 
magnetic impurity remains unscreened for J/D < 0.35. Note 
that Tx\oc changes only weakly with temperature in this case. 
For T <C A it approaches its free value 1/4 (not shown). The 
inset shows the corresponding LDOS where one sees that not 
only the weight is shifted towards the upper half of the band, 
away from the Fermi energy (which is set at quarter filling), 
but also that there is a minimum in the LDOS at the Fermi 
energy, resembling a local pseudogap. 

Free Moment Phase. Next, we use NRG to analyze in detail 
the impurity sites where the renormalization of the magnetic 
susceptibility remains so small that no finite Kondo tempera- 
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FIG. 5: (Color online) The distribution of p c g for W = 2 and var- 
ious lattice sizes. Inset: The magnetic impurity phase diagram in 
the (J, W) plane. The data points correspond to J c = l/p e s miB as 
function of W, obtained numerically for several system sizes. The 
straight line is a guide to the eye. The critical exchange coupling 
expected due to the hard gap caused by Anderson localization, 
[see Eq. $5$] is plotted for comparison (full line). For J > J c , 
all magnetic moments are screened (Kondo phase). For J < , 
all magnetic moments remain free due to Anderson localization, the 
Anderson free moment phase. In the intermediate phase there is a fi- 
nite probability that some magnetic moments remain unscreened due 
to the random occurence of local pseudo gaps. 



ture can be defined. On such sites a magnetic moment remains 
free down to the lowest temperatures, the spin susceptibility 
Xioc diverges and T\\oc{T) approaches a finite value, as seen 
in Fig. [3]3. We employed this criterion to identify free mag- 
netic moments. In Fig. [4] we compare the resulting fraction of 
free moments with that derived from the condition that Eq. (f3]) 
has no finite Tk solution. We see that the fraction of free mo- 
ments obtained from NSE is smaller than the one from NRG. 
This can be attributed to the fact that higher-loop corrections 
tend to lower Tk- Thus, the criterion that the NSE has no so- 
lution gives a lower bound for the fraction of free magnetic 
moments. This criterion corresponds to setting Tk — > in 
Eq. (fJJ. Therefore, free moments exist when the effective, 
weighted LDOS as defined by, 
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is equal or smaller than the inverse exchange coupling 1/J. 
For a fixed realization of the disorder, this yields a lower 
bound for the critical exchange coupling J c below which the 
magnetic impurity remains unscreened. For a clean system 
with a flat band, p c g diverges logarithmically with the num- 
ber of states N or, equivalently, the volume. Then, J c /D ~ 
l/mTV. Since Tk ~ Dexp(—D/J), the condition for free 
moments becomes A > cTk, with a constant c> 1 which 
depends on E-p. In a clean system the level spacing vanishes 
in the thermodynamic limit, therefore there are no free mo- 
ments for any finite J. Can one conclude the same when 
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nonmagnetic disorder is present? No, for the following two 
reasons. First, let us note that for disordered systems with di- 
mensions d < 2, all energy eigenstates are localized in the 
absence of spin-orbit interaction or a strong magnetic field. 
A finite local gap of order A c = l/(v!; d ) develops at the 
Fermi energy, where v is the average density of states and £ 
is the average localization length. Therefore, there are with 
certainty free moments whenever A c ^> Tk, or, equiva- 
lently, J < J C A - D/\nN c , where iV c = D/A c . How- 
ever, in weakly disordered two-dimensional electronic sys- 
tems = g exp(irg) in the absence of an applied magnetic 
field. Here, the dimensionless disorder parameter g = E-pT 
is for quarter filling related to the disorder amplitude W as 
g = 96/ (nW 2 ). Using Eq. © we then find 



J^D [\n(2aW) 2 )+C\ 



(5) 



where C « 0.58. This expression only provides a lower 
bound for J c because both the localization length and the 
LDOS are widely distributed in a disordered system. In Fig. [5] 
we show the distribution of p c ff for moderate disorder and var- 
ious lattice sizes. Remarkably, the point where the distribution 
drops sharply to zero, /9 c ff min , hardly depends on lattice size. 
In the inset of Fig. [5] we plot J c = l//O c ff min as a function of 
disorder for various system sizes, together with J^(W), Eq. 
(01. This plot shows a quantum phase diagram in the (J, W) 
plane with a free magnetic moment phase due to Anderson 
localization for J < J A and a Kondo screened phase for 
J > J C (W). For intermediate couplings, J A < J < J C (W), 
we also find free magnetic moments for all system sizes con- 
sidered. It is well known that a free moment phase exists when 
there is a pseudogap around the Fermi energy [15]. Indeed, as 
seen in Fig. [3}), there is a dip in the LDOS at sites where the 
magnetic moment remains unscreened. In weakly disordered 
metals, the LDOS is correlated over a macroscopic energy in- 
terval of order of the elastic scattering rate 1/r 02611 . Al- 
though these correlations are only of order 1/g and the LDOS 
can fluctuate in energy, they can still cause dips in the LDOS 
within a range of order 1/r. From the numerical results for 
J C (W) (inset of Fig. 5) we infer that the local pseudogaps 
in weakly disordered metals have a power that increases with 
disorder as a ~ W/D ~ l/y 7 ^- 
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